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Recently the ATIC and PAMELA collaborations released their results which show the abundant 
e excess in cosmic rays well above the background, but not for the p. Their data if interpreted as 
the dark matter particles' annihilation imply that the new physics with the dark matter is closely 
related to the lepton sector. In this paper we study the possible connection of the new physics 
responsible for the cosmic e excesses to the neutrino mass generation. We consider a class of 
models and do the detailed numerical calculations. We find that some models can account for the 
ATIC and PAMELA e and p data and at the same time generate the small neutrino masses. 



PACS numbers: 14.60.Pq, 95.35. +d 

A. Introduction 

Recently the ATIC Q and PPB-BETS collabora- 
tions have reported the electron/positron spectrum mea- 
surement up to ~ 1 TeV, with an obvious bump above the 
background at - 300 - 800 GcV and - 500 - 800 GeV, 
respectively. At the same time, the PAMELA collabo- 
ration also released their first cosmic-ray measurements 
on the positron fraction 0] and the p/p ratio [J]. The 
positron fraction shows the significant excesses above 
10 GeV up to ~ 100 GeV, compared with the background 
predicted by the conventional cosmic-rays propagation 
model. This result is consistent with the previous mea- 
surements by HEAT [| and AMS @. Clearly, the ATIC, 
PPB-BETS and PAMELA results indicate the existence 
of a new source of primary electrons and positrons while 
the hadronic processes are suppressed. It is well known 
that the annihilation of the dark matter (DM) can be a 
possible origin for the primary cosmic rays, which can 
account for the ATIC/PPB-BETS and PAMELA data 
simultaneously 0. However, the p/p ratio gives strong 
constraint on the nature of dark matter particles. The 
widely discussed lightest neutralino in supcrsymmetric 
models as dark matter candidate may be difficult to ex- 
plain the PAMELA data due to the constraints @, 0- 
These constraints may imply some special relations be- 
tween the DM sector and leptons. 

Another interesting topic related to the leptons is the 
neutrino mass-generation. In order to understand the 
small but nonzero neutrino masses, required by various 



'Electronic address: bixj@ihcp.ac.cn 
^Electronic address: pgu@ictp.it 

t Electronic address: tli@itp.ac.cn] 

§ Electronic address: xmzhang@ihep.ac.cn 



[MIFP-08-35] 

neutrino oscillation experiments, we can consider the el- 
egant seesaw [l(| extension of the standard model (SM). 
The original seesaw scenario is realized by integrating 
out the heavy particles at tree level [l(J G~H ■ As an al- 
ternative, one may consider the radiative seesaw models 
[12, EH EH EH where the small neutrino masses are gen- 
erated at loop level. 

In this paper we connect the leptonic annihilation 
of the DM to the mass generation of the neutrinos. 
We find that in some radiative seesaw models the neu- 
tral Majorana or Dirac fermions as the DM can dom- 
inantly annihilate into the charged leptons and neu- 
trinos. For proper choice of model parameters and 
boost factor, the leptonic annihilations of the DM 
can explain the electron/positron excesses observed by 
the ATIC/PPB-BETS and PAMELA cosmic-ray exper- 
iments. With small modifications, we might explain the 
INTEGRAL [3 experiment via the exciting DM [l7t | . 
In addition, we propose one kind of models with tradi- 
tional seesaw mechanism (lo| . These models can explain 
the ATIC/PPB-BETS and PAMELA cosmic-ray experi- 
ments similarly, and explain the INTEGRAL experiment 
as well. 



B. Effective Lagrangian 

We start with our discussions with an effective La- 
grangian 



C ^-VtCL,RH X + H.c, (1) 

where i = e,fi,r are the charged leptons, 77 ± is a scalar 
carrying the same charge as £ ± , x denotes a neutral 
fermion. We further assume 77 ± and \ to have no other 
Yukawa couplings with the SM particles. Therefore, in 
case of < m ± , x can dominantly annihilate into the 
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charged leptons since the loop-order annihilations of \ 

into quarks and/or photon are highly suppressed. Fig. 

[1] shows the annihilations of the Majorana x the 
charged leptons while Fig. [2] is for the Dirac X- 



FIG. 1: The annihilations of the Majorana \ into the charged 
leptons through the exchange of r] . 



The cross sections of x annihilating into the charged 
leptons are given by 
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where x is a Majorana or Dirac particle, respectively. ing particles in their cms system. Up to 0(v 2 ), the cross 
Here v is the relative velocity between the two annihilat- sections (f2a| and (|2bj) can be respectively simplified as 
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with the definition, 
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By inputting v ~ 10 3 (the average velocity in our 
galaxy) and r = 0, we further derive 



(4) 
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(? a p v ) < 5.9 x 10~ 26 cm 3 sec" 1 ( 700 GcV J \y a \ 2 \ yi3 \ 2 inDiraccase, (5b) 
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FIG. 2: The annihilations of the Dirac \ m t° the charged 
leptons through the exchange of rj . 

where y e .n,T should be smaller than v^Li" for a perturba- 
tive theory. 

In the following, we will embed the effective La- 
grangian Eq. ([T]) into some radiative seesaw models where 
the tiny neutrino masses arise from loop contributions. 
And we can embed the effective Lagrangian Eq. (Q]) into 
the models with traditional seesaw mechanism as well. 
For our convention, we denote the SM lepton doublets, 
right-handed leptons and Higgs doublet as ip L (2, —h) = 

{v La ,ll a ) T , e R Jl,-l) and <t>{2,-\) = (0°,0-) T , re- 
spectively, where their SU(2)l x U(1)y quantum num- 
bers are given as well. 

C. Model I with Majorana Neutrino Masses and 
Mixings from Radiative Corrections 

We consider a Majorana radiative seesaw model [TJ|, 
fl3j where one introduces three right-handed neutrinos 
X R . and a SM scalar doublet r\ with SU(2) L x U{1) Y 
quantum numbers (2, —5). In this model, Xr, and V 
are odd under a Z 2 symmetry, while the SM fields carry 
even parity. Respect to this Z 2 , the scalar doublet r\ will 
not develop nonzero vacuum expectation value (VEV). 
Therefore, as shown in Fig. [3J the neutrino masses can 



FIG. 3: The loop diagram for generating small masses of Ma- 
jorana neutrinos. 

be radiatively generated through the following Yukawa 
couplings, the Majorana mass term of the right-handed 
neutrinos, and the quartic interaction of scalars 

-CD ^m x .x%'.x Ri + Vaitp^VXRi + A(0 f ??) 2 + H.c. .(6) 

Here we have conveniently chosen the base, in which the 
Majorana mass matrix of the right-handed neutrinos is 
diagonal and real, i.e. m = Diagjm , to Vo ,to }. Ob- 
viously, the above Lagrangian in Eq. (|6J) can induce the 
effective theory in Eq. JT|). If the lightest \ i s lighter 
than r\ R = -^Rc(r/°), j] I = -^Im(?7°) and rf 1 , it will 
have no decay modes and will dominantly annihilate into 
the charged leptons and the neutrinos. The total cross 
section of the lightest x m to the charged leptons and 
neutrinos is given by 

(av) < 3.0 x lO^Wsec" 1 (^^) 2 [(^kf -(7) 

Here m ~ m ~ m 2 + has been adopted since we are 

interested in the case where the Yukawa couplings y ai are 
large enough for a desired cross section and the quartic 
coupling A is small enough to guarantee the smallness of 
neutrino masses. 

Because the two light right-handed neutrinos xi and Y2 
can be degenerate, we might explain the INTEGRAL [16[ 
experiment via the exciting DM [TjJ. To explain the 
INTEGRAL experiment by producing X2 via the chan- 
ne l Yi Yi — > X1X2 through the tree-level and ladder dia- 
grams [13 , we introduce a SM singlet S that is Z 2 even. 
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The discussions are similar to the following subsection 
D except that wc will have additional terms in the La- 
grangian in addition to that in Eq. (fT4|) 



H.c. 



(8) 



where Ai are mass dimension-one parameters. In short, 
this approach needs some fine-tuning. 

To avoid the fine-tuning, we can consider the large dis- 
crete symmetry T in the dark matter sector, for example 
Z4 in Ref. [17| . or as discussed in subsection D. Note 
that the Majorana masses for the right-handed neutrinos 
are forbidden by T since it is not Z 2 . We will need two 
SM singlet Higgs fields: one generates the large Majo- 
rana masses, while the other generates the small mass 
splitting around 1 MeV and has Yukawa coupling about 
0.18 or larger to the DM and exciting DM fields so that 



we can produce \i [13 • Because this kind of model is 



similar to the models in the following subsection C, we 
will not study it in detail here. 



D. Model II with Dirac Neutrino Masses and 
Mixings from Radiative Corrections 

Let us consider a Dirac radiative seesaw model [l4| by 
extending the SM with a SM scalar doublet 77(2, —3), a 
complex scalar £(1,0), a real scalar cr(l,0), three right- 
handed neutrinos i/ R (1,0) and three Dirac fermions 
XLi^i 1 ^)- In tnis model, there is a U(l) 

d gauge sym- 
metry under which £, v R and x R carry the quantum 
number 1. In addition, 77, a and Xl r are °d u while 
all of other fields are even under a Z 2 symmetry. More- 
over, i/ R and Xl- _r ; are arranged for the same lepton 
number with the SM leptons. Conserving the U(1) D , 
Z 2 and lepton number, we can write down the following 
interactions 



£ y a ii ) L a , nx c L i + h a i(Ti' Ra x Ri 

J rfiiCX R ^XL j + /W</> + H.c. 



(9) 



Since 77 and a are forbidden to develop the nonzero 
VEVs as a result of the Z 2 protection, the neutrinos can 
only obtain their small Dirac masses via the one-loop dia- 
gram, as shown in Fig. |4j It is straightforward to see the 
effective Lagrangian in Eq. (fT]) can be embedded into this 
Dirac radiative seesaw model. Note there exists a mixing 
between £ and <j>° , which will lead to the co-annihilation 
of the dark matter into the quarks. But (£) should be 
about one order larger than (</>°) for giving the dark mat- 
ter a mass close to the TeV scale. So, the £, — (j>° mixing, 
of the order of 0((</>°)/(£)) ~ 0.1, could have no signif- 
icant implications. For convenience and without loss of 
generality, we, in the following, shall choose the basis in 
which the Dirac mass matrix of the singlet fermions is di- 



Similar to that in the Majorana radiative seesaw model, 
here the lightest x can also dominantly annihilate into 
the charged leptons and the neutrinos if it is lighter than 
r]j = -^Im(?7 ), ri^ and the other two physical states 

p 1 2 (combinations of = -i=Rc(7/ ) and a). The total 
cross section of this lightest x m to the charged leptons 
and neutrinos is given by 



(av) < 5.9 x 10 



" 26 cm 3 sec" 1 



700 GeV 



{[Tr(/,t/ l ) rt ] 2 + 2[Tr(7 / t 2/ y 2 } . (10) 



has been adopted by 



Here m 2 Pi ~ m 2 2 ~ m 2 ? ~ m 2 j 
taking small trilincar coupling /1. This is consistent with 
our scope for large Yukawa couplings y a i and h a i and 
then a desired cross section. 

Naively, we may guess that wc do not need to introduce 
another SM singlet Higgs field £' to explain the INTE- 
GRAL experiment. However, the Dirac masses for Xlr 
are generated from the terms /yCxjjXz,, and we need 
two degenerate Dirac fermions Xlr m ^ nc ma ss basis. 
Because the masses for dark matter fields is about 600- 
800 GeV and the mass splitting for the two light Dirac 
fermions is about a few MeVs, we obtain that the Yukawa 
coupling constants for £x'lXr and £,XlXr are about 10~ 5 
unless we fine-tune the coefficients fy. Thus, it is dif- 
ficult to generate x'lr through the tree-level and lad- 
der diagrams, and then explain the INTEGRAL exper- 
iment [131 ■ Therefore, we introduce another SM singlet 
Higgs field £' with VEV about 10 MeV. f has U(1) D 
charge +1 and is Z2 even. The additional Lagrangian is 



C 



+f(^(^) + ^WO 



Ai 



^OW^+H.c, 



(11) 



agonal and real, i.e. m x = /(£) = Diag{m xi , 



To explain the INTEGRAL experiment, we need to 
be about 0.18 or larger (l7j . 



E. Model III with Neutrino Masses and Mixings 
from Traditional Seesaw Mechanism 



We consider the generalized Standard Model with two 
or three heavy right-handed neutrinos. The neutrino 
masses and mixings are generated via the traditional see- 
saw mechanism, and the observed baryon asymmetry is 
explained via leptogensis. 
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FIG. 4: The loop diagram for generating small masses of Dirac 
neutrinos. 



First, we consider the model that can only explain the 
ATIC and PAMELA experiments. We introduce dark 
matter sector that has a SM singlet scalar field E with 
SU(2)l x U(1)y quantum numbers (1,-1) and a Majo- 
rana fermion %. To have a stable dark matter candidate 
X, we consider the Z2 symmetry under which only E and 
X are odd. The renormalizable Lagrangian for the dark 
sector and its interaction with the SM is 



1 _ 1 2 ~t~ 



Ai 



{EflEf 



+^{&E){<t?<t>) + (tf e e i R Ex + K.c) ,(12) 

where m x and rriE are mass terms, and Ai and y\ are 
Yukawa couplings. Also, Ai and A2 are real, while y\ can 
be complex and then violate CP symmetry. To have \ 
as a dark matter, we require that m x <mj. 

Second, we can generalize the dark matter sector in 
the above model to explain the INTEGRAL experiment. 
In the dark matter sector, in addition to E as above, we 
introduce a SM singlet scalar field S, and two Majorana 
particles \i an d X2 respectively with left and right chiral- 
ities [p?! . And the discrere symmetry in the dark matter 
sector is Z4. Under Z4, the fields in the dark matter 
sector transform as follows 



Xi,2 - e^ 2 Xl ,2 , E - e-™ / 2 E , S 
The renormalizable Lagrangian is 
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m D XiX2 

^2 rcl \ n2 
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(14) 



In the basis Xi X* 
for x's is 



l/y/2(xi T X2), the mass matrix 



M = 



y+S-m D y-S 
yS y + S + m D 



(15) 



where y± = \{yi ± 2/2)- Thus, at the leading order, the 
Dirac fermion can be decomposed as two degenerate Ma- 
jorana fcrmions x an d X* ■ If the Z4 symmetry is broken 
weakly to Z2 by giving S a VEV, we expect these states 
to be split by a small amount 5 = 2y + (S). Suppose 



< 0, we obtain that (S 



and the mass 

splitting is then just 5 = y+\ms\/y/X2 To explain 

the INTEGRAL experiment, the mass splitting is about 
a few MEV. Note that the couplings y± can be as small 
as around 0.18, and then the VEV of S should be around 
10 MeV. Thus, the mixing between the Higgs field <fi° and 
S is very small from A5 term, and then the anti-protons 
from the XX annihilation arc highly suppressed. 

Interestingly, this kind of models can be naturally cm- 
bedded into the supersymmetric flipped SU(b) x U(i)x 
models [l|[ and the string or F-theory derived flipped 
SU(5) x U(l)x models [19(. An important difference of 
this kind of models from the other two models is that it 
predicts no neutrino signals from DAI annihilation. The 
models I and II give large neutrino fluxes as well as the 
positron flux. Detecting neutrino fluxes from the Galac- 
tic Center may discriminate these models. 



F. Phenomenology 

In this part we study the dark matter models given 
above in detail to explain the ATIC and PAMELA data. 
The propagation of charged particles in the Galaxy is 
calculated numerically in order to compare with data 
measured at the Earth. The charged particles propa- 
gate diffusively in the Galaxy due to the scattering with 
random magnetic field [2p| . The interactions with inter- 
stellar medium will lead to energy losses of the primary 
electrons and positrons. In addition, the overall convec- 
tion driven by the Galactic wind and reacceleration due 
to the interstellar shock will also affect the distribution 
and spectrum of electrons. The propagation equation can 
be written as \2n 



dip 
~dt 



Q(r,p) + V-(D xx Vi>-V c ip) + 



d_ 

dp 



d_ r 

op 



(16) 



where ip is the number density of cosmic ray particles per 
unit momentum interval, Q(r,p) is the source term, dis- 
cribing the primary particles injected into the interstellar 
meadium, D xx is the spatial diffusion coefficient, V c is 
the convection velocity, D pp is the diffusion coefficient 
in momentum space used to describe the reacceleration 
process, p = dp/dt is the momentum loss rate, mainly 
induced by synchrotron radiation and inverse Compon 



6 



scattering. In the work we solve Eg. (|16|) numerically 
adopting the GALPROP package [22[ ■ 

For DM annihilation, the source term of electrons and 
positrons have the form 



Qa(t,E) =BF 



< av > A p 2 {r) dN(E) 



dE 



(17) 



where, BF is the boost factor, < av >a is the DM anni- 
hilation cross section, ttidm is the mass of DM particle, 



and 



dN(E) 
dE 



is the electron/positron spectrum from one 
pair of DM annihilation, p(r) is DM density distribution 
in the Galaxy. In this work, we adopt the NFW pro- 
file for DM distribution with the local DM density 
Po = 0.3GeV/cm 3 . 



ATIC f- 
H'B b: b I >• 

HESS - 
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PAMELA i— •— i 
HEAT9495 

HEATOO i — i — i 

background 

ba ckg round +62 OGeV(aV) . 

:aci:onuuriu+780GeV(e,mu,taj) 







FIG. 5: Top: Electron/positron spectrum including contribu- 
tion from DM annihilation compared with the ATIC/PPB- 
BETS data; Bottom: e + /(e~ + e + ) including contribution 
from DM annihilation as function of energy compared with 
the PAMELA and HEAT data. Two models are considered: 
in one model DM mass is 620 GeV with e + e~ being the main 
annihilation channel, while in the other model DM mass is 
780 GeV with equal branching ratio into e + e~, p + p~ and 

T + T~ . 



780 GeV and assume DM annihilates into e + e~, p + p~ 
and t + t~ with equal branching ratios. We can see 
that both models can give good fit to the data. The 
value BF < av > is taken as 7.2 x 10 _24 cm 3 s _1 and 
1.8 x 10 _23 cm 3 s _1 respectively for the two models, cor- 
responding to the boost factor 240 and 600 respectively 
if we assume DM produced thermally with natural value 
< <tv >= 3 x 10 _26 cm 3 s _1 . In order to explain the 
ATIC/PPB-BETS and PAMELA/HEAT data we con- 
sider the scenario where the DM is produced by some 
nonthermal mechanism (24j. According to Eqs. (|5a[) and 
(|5b[) with the Yukawa couplings y a ^ < V^Lr, we can see 
that: (i) for the Majorana fermion as dark matter, the re- 
quired boost factor should be larger than O(10 7 ), which 
is in lack of reasonable explanation; (ii) for the Dirac 
fermion as dark matter, the boost factor can be absent 
or can be as small as 0(10), which may be due to the 
clumps of DM distribution [25[ . 



G. Summary 

The recent observations by ATIC/PPB-BETS and 
PAMELA imply that DM should dominantly annihilate 
into the leptons. This phenomenon inspires us to con- 
nect the leptonic annihilation of the DM with other lep- 
tonic processes, such as the neutrino mass-generation. In 
some radiative seesaw models, where the neutrinos ob- 
tain their small Majorana or Dirac masses at one-loop 
order, DM can naturally have dominant channel anni- 
hilating into the leptons. We show that the Dirac ra- 
diative seesaw model can account for the observed elec- 
tron/positron excesses even if we don't introduce large 
boost factor. However, the Majorana case is strongly 
disfavored due to a unreasonably large boost factor to 
account for the PAMELA/ ATIC data. With small mod- 
ifications, we might explain the INTEGRAL experiment 
via the exciting DM. Furthermore, we constructed one 
kind of models with traditional seesaw mechanism. These 
models can not only explain the ATIC/PPB-BETS and 
PAMELA cosmic-ray experiments similarly, but also ex- 
plain the INTEGRAL experiment. 
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In Fig. [5] we show our results compared with the 
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